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Annihilation process of proton and antiproton to guark antiquark pair accompanied by emission 
of two additional gluon jets with intermediate state of vector meson is considered. Strong azimuthal 
correlation is revealed of two gluonic jets, effectively created in the same plane. Some applications 
to cosmic ray events and to LHC experiments are discussed. 



^ , I. INTRODUCTION 

' Now is widely believed that hard processes with hadron interaction at high energies can be described in frames 
" of QCD. Contrary to photons in QED the QCD gluons can interact between themselves. The experimental check 
of nature of gluons become an important problem up to now. The dominant role of branching gluons was utilized 
^ in construction of the gluon dominance model (GDM) describing the events with large multiplicity in high energy 
collisions of leptons, (anti)protons and ions [l]. Another indirect indication can be found in measuring the azimuthal 
^\ ' correlations between gluon jets emitted by color quarks created in some hard hadronic process. It is the motivation of 
this paper to show that such kind of correlation take place in particular in process of annihilation of the high energy 
(-H proton and antiproton to quark-anti quark with emission of two gluons. 

O Years ago in the paper of one of authors |2] a problem of inelastic quark form factor was considered. The con- 

I . sideration was performed in the so called double-logarithmical approximation (DL) when the leading terms of order 
{as /n)^ {\n{Q^ /m^))'^ was taken into account (here — —q^ is the square of 4-momentum of virtual vector particle, 
(— I , m is the effective mass of the gluon jet and as is QCD coupling constant). The one- loop radiative correction to the 
I— I ■ inelastic form- factor with a single hard gluon emission was found as well as the contribution from the channel of two 
hard gluons creation. It was confirmed, in particular, the S. Adler statement about cancelation of all DL enhance- 
ments when considering two-loop virtual corrections together with contribution of inelastic channels. Compared with 
QED case apart from some complication connected with the color properties of quarks and gluons, some additional 

■ kinematical regions arising from the existence of three gluon interaction vertex become to be important. It is the 
QQ ' motivation of this paper to generalize the results obtained in 1978 year to the annihilation channel. Another reason 
T— I . is to search the possible relation with some azimuthal correlation of particles observed at LHC [s*] . 

■ Some additional speculation can be done, a prediction about the character of energy distribution of cosmic rays of 
. high energy Really, taking into account the strong azimuthal correlation of jets created in high energy proton 

T-H ■ with nuclei in atmosphere (all the jets are in the same plane) a specific distribution of the spots of sets of the 
^-H secondary particles of high energy can be measured: all these spots must be located along some lines-intersection line 
of production plane and the surface of the Earth. 

Process of annihilation of proton and anti-proton to a vector meson with the subsequent creation of quark-anti-quark 
pair and two hard gluon jets 
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p{p+)+P{P-) ^ v{q) q{q^) + q{q+) + g{ki) + g{k2) (1) 

in Born approximated is described by a set of Feynman diagrams of two kinds. One of them corresponds to emission 
of both gluons by quark and anti-quark and another contain the specific for QCD three gluon vertex, when two gluons 
arise from a single gluon, emitted by quark or anti-quark. We will use Sudakov parametrization for 4-momenta of the 
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problem, introducing two light- like vectors 

h = piPi +P2ai + kij_,ki±pi^2 = 0; kfj^ = -fcf < 0; 

2 2 / ^2 1-/3 1-/3 

Pi=P2 = 0;s= {q+ +q-) ,pi = q- - y^I+'^P^ = 1+ - YT^'^^' 

ql = ql=M^='-{l-n 

k1 = sUiPi - I? = . (2) 

Here M is quark mass and m is the mass of vector particle - gluon (or the invariant mass of the relevant giuon jet). 
In terms of Sudakov variables the phase volume of two real soft gluons can be written in form 



(Pki(Pk2 _ ,4; xn2 2\j4i 2 2\ 
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(fki5{ki - m^)(fk25{k^ - m^) 



-daid^i(f kiS{sai/3i — /c^ — m'^)-^da2dl32d? k25{sa2p2 — — m^) = 

(i)^^daidMMi - l)^da2dM{a2f32 - 7), (3) 

with 7 = mVs << 1. For the ratio of matrix elements squared and summed on final state quantum numbers was 
obtained in 

~ ^^^1 ^ ^n^\l3i + P2){ai + a2y 2N ^ 

1 1 1 1 

4iV2 ^^^^2^ ai^2 (^1 +;32)(ai +02) ^ aia2/3i/32 

H ^ i \ \ 1-4+ — + — + — + —III 

4((fci+A:2)2-m2)V2/3i ai/32 (/?! + /32)(ai + "2) ^ ai /32 /3i 

= ^+77 A-^ 9 5- (4) 

(fci + fc2)2 -m2 ^ ^ 

A'^ = 3 is the rank of color group SU{N). Main contribution arises from two kinematical configurations in 4-dimensional 
manyfold R: j < ai, l3i < l,ail3i > 7 with the additional conditions in the region a) 

ai»aj, I3i»l3j, i,j = l,2 

and region b) 

1 >> |ai/32 - a2l3i\ » 7. 
Contribution of first kind kinematics, region a), is 

^^U, = |(^)V4. . = c, = ^ (5) 

where as is gluon coupling constant and dao is the differential cross section of annihilation of proton and anti-proton 
to the pair of quark and anti-quark (see Appendix A). 

But it is not a total answer. The double- logarithmic type contribution arises as well from the region b). Really the 
denominator of the virtual gluon propagator 

(fci -I- /c2)^ — — 2kik2 +m'^ — s(ai/32 + a2/3i) + — 2|fci||fc2| coscj) = s{a ~ 6cos0), 

k\ = sai/3j - m^. (6) 

Being averaged on azimuthal angle it has a form 



277 a — 6 cos ^ \fa}~^\p- 
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Using the on mass shell conditions for gluons we have 

1 



(8) 



Va^ - 62 - a2/3i)2 + 47[ai/?2 + ai/3i + a2/32] 

We see that in region b) it is presented a quite different region of DL contributions to the cross section. 
Performing the integration it was obtained ( [H expression (17)): 

Color structure containing cy = N violates Poisson distribution in emission of two colored vector particles 

= _a"e^ a^— . (10 

dcTo n\ ZTT 

Here factor with exp{—a) takes into account the contributions from virtual vector meson emission. 
We suggest to measure the azimuthal correlation in form 

rf^(2)(0) 



AW 



J da(^){(j))d(j) 



4 + ScFCvTT^^ 4 + 72n^^ 



Function Z{p, (j)) satisfies the normalization condition J A{4))d(f) = 1 
To obtain Z{p, cj)) we simplify the expression a ~ bcos(f) in such a way: 



(11) 



a — b cos (j) — b^ -\- a^c/)^ 

Performing the a2 integration we obtain 

da2 47rai/32 1 



1 2a 

a = 2ai/32. (12) 



a-bcos(j) I3{ y/R 
The quantity R have a different form depending on integration regions. We obtain (see Figs. 2-4): 



(13) 



1 J ai J V-Ri J VR2 J Oil J Pi J yjRi 

/3l/ai 



'^^X14) 



with 



4to2 q,^ ai/320^2 

-n-i = -5- + 4(— ) ; 

s Pi Pi 



We use the integrals 



j d02 I3i ln(l//3i) 



R2 2ai 



^1 ^"i/^i 

/3l/ai 

d/32 /3i ln(l/ai) 



/Ih 2ai 

SQl/3l 



(16) 
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and 



dh Pi 

X 



/i?r 2ai|0| 



{ln(i|0| + v/T+iV) - ln(|0| + + i = 7^1^- (17) 
To express the result in terms of one-fold integrals we use 
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\/7 7//3l 

1 Pi 

^ln(l//3i) / '^F{t)^l f dy{p-yfF{eMym)\ 
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(18) 



\/7 7//3l 



The result is 



p 

Z{p, 0) = J /(P - 2/)' . / + 

P 

l-J{p-y) ln(|0|e^/2 + yTT0^)dy + 0(p2). (19) 




Function Z{p, </>) has a delta-function type behavior, concentrated in \(p\ ^ ^ s/m? with Z(0) « Sy's/m^ and besides 
/ Z(</.)# = (l/24)p^ Here ^/^^ « e « 1. 

— e 

This function is presented in Figure 1 for p = 10. 

II. DISCUSSION 

Let remind the three dimensional picture in center of mass of initial particles. Hard quark and anti-quark of final 
state move back to back at large angles and emit two gluons. One of them moves close to quark direction another - 
close to anti-quark ones. Differential cross section on the azimuthal angle between the planes containing quark and 
the relevant gluon and the plane containing anti-quark and another gluon will have a isotropic part and one with 

sharp dependence distributed close to the value \(t)\ « 0{\J~^). 

Energies of gluons are small compared with energies of quark-and anti-quark. This statement follows from the 
fact that the main,(^ p^) contribution follows from the 4- vectors of gluons polarization in kinematics of isotropic 

contributions to the cross section are situated in the plane of quark and anti-quark 4-momenta ei ~ — ^ ■f^P2 + \J '^Pi ; 

whereas the ones, responsible for azimuthal correlation are essentially situated in the plane transversal to quark 
momenta. For the case of production of additional soft quark-anti-quark instead of two gluon the DL enhancement 
factor do not appears. 

The DL enhancement phenomena will take place for gluons emitted by any pair of colored fermions which have 
large invariant mass ^/s » m. However the explicit form of (^-independent and (/)-dependent parts of matrix element 
square (see (4)) will depend on mechanism of quark creation. In particular one can consider the peripheral mechanism 
of quark-anti-quark creation in collisions of hadrons by two reggeized gluons. 

From the comparisons with experimental topological cross sections of pp interactions and pp annihilation at the 
energy region close to U-70 accelerator (IHEP, Protvino) the estimations of the MGD-parameters were obtained [l|. 
The mean number of hadrons formed from one gluon source at hadronization stage is added from the charged and 



5 

A((/j) 




-300 -200 -100 100 200 300 

Fig. 1: the azimuthal distribution function A((j!>) for j s — 10~*, — rj^/mFJs is presented. 
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Fig. 2; Integration domain in qi, /?i plane. 

neutral components (pions predominate): 7Ttot=l. 63+1. 01=2. 64. The invariant mass of such gluon sources may be 
determined as m— 2.64 x 0.139= 0.37 (GeV). At accelerator U-70 energy GeV. The parameter p will be 

equal to, p = In -^=6.9. At that p the noticeable peak may appear in the angular dependence too. It looks also 
as 1-GeV g-jet at low threshold of ISR-energy, y^=32 GeV. To describe the topological cross section widening at 
ISR energies the gluon fission was included at QCD-cascade stage in GDM py. One can make next assumptions. To 
reveal the angular distribution peak appeared from gluon fission at the energy lower than LHC region it is enough to 
choose g-jets with smaller invariant mass m or determine mass of g-jet by means the peak. Also the estimation of the 
invariant gluon mass is compared to the mean transverse momentum of secondary particles at the relevant energies 
by the surprising way and this can be clue to ridge phenomenon understanding. 

III. APPENDIX A 

Matrix element of process p{p+) + p{p-) v(q) — > q{q-) + <jiq+) have a form 

^ 9Vvv9v,,F,y) -^^_^^^^^^^^^^ ^ «(g-)7.i'(9+) (20) 

q IVly 

with phenomenal approach for the quark form factor and gvpp = '^gvqq ~ 9, is dominantly w-meson. 

Differential cross section have a form (c = cos 9, 9 is the center of mass angle between the directions of motion of 
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Fig. 3: Integration domain for I region (see Fig. 2 )in 02, /32 plane. 




Fig. 4: Integration domain in 02, /32 plane for II region (see Fig. 2). 



initial proton and the positively charged quark with momentum q^) 

dao 1 



27rs/3(s - M2 )2 

i?o = Ri\F^? + T^[\F^? - \Fe?)R2, (21) 

1 — T 

with iV = 3 is the number of quark colors 

F„ = Fi + F2; Fe = Fi + TF2, r = ^ (22) 

and 

i?i = i2 ^ tt2 + 4,5(m2 + ml) - 2{M^ + mlf; 

R2 = tu + s{ml ~ {M^ + ml f . (23) 

Here we use the kinematic invariants 

s^{p++ P-f\ t={p+- q+f- u+{p+- q-f; s + t + u = 2{M^ + m^), (24) 

M, m^-proton and quark masses. In paper Q was argued that dominant contribution is provided by l^-omega meson, 
and, besides F^ = 0. In paper the reasonable model for Dirac formfactor of proton in annihilation channel was 
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suggested 
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